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Abstract 

We study a two-dimensional bosonic field theory with a random defect line. The 
theory has a background field coupled to the field variables at the defect line, which 
renders the model non-integrable. However, as the background field is random, 
and the disorder is implemented through the replica trick, the model becomes inte- 
grable, allowing us to use the form-factor method to compute the exact correlation 
functions of the quenched model. 
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1 Introduction 

Integrable field theories in two-dimensions [1] are at the same time useful from the the- 
oretical poin of view as well as to applications to statistical mechanics and condensed 
matter problems [2, 3]. Recently there has been a considerable amount of work on two- 
dimensional quantum field theories with boundaries and defect lines [4]- [16]. The interest 
on this type of theories is due to the fact that they can model such systems as quantum 
impurities, for example. It has been shown that if one wants the theory to be integrable 
even after the defect has been introduced, then the bulk theory has, necessarily, a con- 
stant S-matrix, being basically that of a free fermion or free boson [6, 11]. This restricts 
quite severely the type of interactions one can use at the defect line in order to mantain 
integrability. 

In this paper we study a two-dimensional free boson with a disordered defect line. 
Before taking the disorder into account the theory is non-integrable, but after disordering 
the line, through the replica trick [20] , we obtain an integrable field theory for all values 
of n, the number of replicas. Another interesting phenomenon that happens concerns the 
stability of the vacuum of the theory. As it has been shown in [6] there is a minimun 
value for the coupling constant, below which the theory becomes unstable. We will see 
that this happens for every finite value of n, but that in the limit n — * the theory 
becomes well-defined for all values of the coupling constant (which is the width of the 
disorder A). 

This paper is outlined as follows. In section 2 we introduce the model to be studied 
and discuss some of its properties. In section 3 we apply the replica trick and obtain the 
effective field theory with finite n. The reflection-transmission algebra is briefly reviewed 
in section 4, and in section 5 we compute the reflection and transmission amplitudes. We 
show how to compute correlation functions using these amplitudes in the next section, 
and comment on the analytical structure of the theory as n — > 0. Our conclusions are 
presented in section 7. In the appendix we collect some of the useful formulas to compute 
correlation functions for the disordered theory in the n — > limit. 

2 The Disordered Bosonic Model 

Consider a two-dimensional free-boson Euclidean quantum field theory (qft) with a line 
of defect, given by the following action 



where £<f(0, <90) is the defect interaction lagrangian. It is clear that, even though the 
model is obviously solvable in the bulk, we have to choose judiciously the defect interac- 
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(2.1) 



tion in order to mantain integrability. In principle we could have chosen different qft's 
for x < and x > 0, as have Corrigan et al done in classical context [14, 15], but in this 
case we will consider the same theory on both sides of the defect. 

If we couple this theory to an external field with a defect action dt h(t)(f)(0,t), 
then the theory is clearly non-integrable, since the external field can absorb or produce 
new particles. We spoil the fact that in an integrable model there is no particle production 
in a scattering process. 

The equations of motion for this model are easily derived to be the usual Klein-Gordon 
equation, plus a source term 

d 2 <P + m 2 <p = 8(x)h(t) . (2.2) 

The coupling with the external field is, therefore, taken into account through the bound- 
ary condition. In order to consider the (quenched) disordered theory we are going to use 
the replicated action, as discussed in the next section. 



3 The Disordered Action: Replica Trick 

When considering quenched theories we need to compute the average of the Helmholtz 
free energy, T . The technical complication that arises immediately is how to do that, since 
we need to compute, essentially, the average of the logarithm of the partition function 



Z. This may be acomplished through the replica trick: we write the average InZ = 
lim n ^o z "~ 1 . By taking the average before the limit n — > 0, we arrive at the replicated 
action. Some of the formulas showing how to compute correlation functions from this 
replicated action are shown in the appendix. 

We consider a gaussian distribution for the background field h. This implies that 
(h(t)h(t')) = A5(t — t') and (h(t)) = 0. All the higher moments can be calculated from 
these two expectation values. 

The replicated action will be given, then, by 

Z$ = fflWkf Vhexp (- £(I / d 2 x \{d<p k f - m 2 ^} + f°° dt h{t)<t> k )- 

k=l \ k=l 1 J J - co 

-klZ ■ (3 - 3) 

Performing the integration over the background field, we obtain 

4? = / ft %exp E( / d 2 x [{d<p k f - m 2 ^]) + | £j dt (± 



k=l \ " k=l " " k=l 



(3.4) 
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This partition function defines the action we are going to study in the next sections. We 
can readily see that the model has become integrable: the defect interaction becomes 
quadratic in the replica copies. An analogous mechanism takes place in the case of the 
ising model with a random boundary magnetic field [17]. 

4 Reflect ion- Transmission Algebra 

Similarly to the case of two-dimensional boundary integrable field theories [4], we can 
start from a given S'-matrix, and look for the constraints imposed on the reflection and 
transmission amplitudes. In what follows we will consider diagonal S'-matrices, for sim- 
plicity. We also assume that the defect has no internal degrees of freedom, the inclusion 
of which is rather elementary. Let us denote the reflection and transmission amplitudes 
by Rij(0) and Tij(6), and let Zi(8) be the Faddeev-Zamolodchikov operator that creates 
an asymptotic state for particle i with rapidity 9. 




Fig 1. Graphical representation of reflection and transmission amplitudes. 

Since the model we are considering does not break parity, we do not have to make a 
distinction between amplitudes of particles reflecting or scattering from the right or left 
side of the defect. Let us introduce the operator T> for the defect. The operators Zi{9) 
and T> obey the following algebra 

Zi(0)V = Rij(0)Zj(-0)V + Tij(0)VZj(9) , 

VZi{9) = Rij^VZji-e) + T lj (-9)Z J (9)V . (4.5) 
Applying this algebra twice in VZi(9), implies the following unitary conditions 

Ri k (9)R kj (-9) +T ik (9)T kj (-9) = 8{ , 

R ik (9)T kj (-9) + T ik (9)R kj (-9) = . (4.6) 

Since we are considering diagonal S'-matrices, the Yang-Baxter equation is trivially sat- 
isfied. 
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5 Reflection and Transmission Amplitudes 

Once we have the equations of motion of the disordered model (3.4) we can use a mode 
decomposition of the left and right fields in order to solve them and find the reflection 
and transmission amplitudes. 

The k-th copy of the field in the replicated action, 4>^ k \ can be written as a sum of 
fields to the left and to the right of the defect 

4 k) (x, t) = 9{x)(f)f (x, t) + 9(-x)<f) {k) (x, t) , (5.7) 

where 9(x) is the Heaviside function. 

Let us introduce the mode decomposition for each replica as follows 

(/) { ±\x,t) = [ + °° — ( a { ±\9) exp(-im(t coshO - xsmhO)) + c.c.) , (5.8) 

where 6 is the rapidity variable, meaning that energy and momentum are written as 
( e >p) — (racosh#, msinh#), and the coefficients a±\0) and a±(9) are the creation and 
annihilation operators, respectively, and they satisfy the usual commutation relations 

[ai\e 1 ),a^ > \e 2 )} = 2*6*6(6! - 9 2 ) , (5.9) 

and all the remaining commutators vanish. 

Substituting the decomposition (5.7) into the equation of motion, we obtain the fol- 
lowing boundary conditions, by integrating around x — 0, and using appropriate regu- 
larizations for the Heaviside and delta functions 

A n 

^(^(0,0-0^(0,0) = 7 E(#(M) + 0^(0,0) , 

4 3=1 

0^(0,0 =0-W(O,O ■ (5.10) 

We see then that the different replica copies interact through the defect only. These 
boundary conditions among the field components become conditions among the modes 
in (5.8) 

• A n 

[af\B) - af\-6) - a^(-e) + a^(-9)) = ~ 4m ^ (g) E(ag )f (9) + a™ (-6) + 

af \9) + af\-9) = a^\9) + a (k) \-9) , (5.11) 

and all the remaining commutators vanish. It is easy to see that these equations agree 
with the ones in [6] if we take n — 1 and make the correspondence A = g, where g is the 
coupling constant of the defect for the bosonic theory studied in their paper. 
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These boundary conditions can be summarised into one matrix equation as follows 

M + [6]V[6] = M-[6]V[-9] , (5.12) 
.. a^\9) a+^(0)\ is a 2n component vector, 



where V[9} T = [a ( l )j (9) aj )f (-9) 
and the matrix M± [9] is given by 



M±[6] 



±1 










±1 












1 ± h{9) 1 ± h{9) 
±h{9) ±h(9) 



±h{9) ±h{9) ±h(9) 
l±h(9) l±h(9) ±h(9) 



(5.13) 



where the first n rows are made out of ±l's and the last n rows of 1 ± h and 
h{9) = iA/4m sinh We can find the reflection and transmission amplitudes 
the algebraic method. First, recall that we have a linear relationship among the 
and annihilation operators that act on the right or left of the defect, as follows 



±h, and 
by using 
creation 



a 



(i)t 
+ 



(0) 

(-0) 



aL n)t (0) 



-0)1 



R\{9) T}(9) 
Tl{9) R\{9) 



T 2 n {9) R 2 n {9) 



T l n {9) 
Ri{0) 



R 2 n(0) 

T 2 n {9) 



Rl(0)\ 



t (1)t ( 



(0) 



(0) \ 



(5.14) 



The physical interpretation of the Rl(9) and T-(9) is quite simple. The Rl(9) are the 
amplitude of probability that an incoming particle of type % moving to the right (left) 
with rapidity 9 (— #)reflects to the left (right) at the defect and becomes a particle of type 
j. Similarly, the T-(9) represent the amplitude of probability that an incoming particle 
of type i moving to the right (left) with rapidity 9 (—9)is transmitted to the left (right) 
of the defect and becomes a particle of type j. 

Using equations (5.12) we see that the reflection and transmission amplitudes are 
obtained from the matrix M±[9] by computing R[9] = M+^iWl [6>]. 

Using (5.13) we can solve for the Rij(9) and Tij(9). This is most simply done by 
noting that the inverse of M+ l {9) has a structure similar to M+(9) 



M-\9) 





1 







-1 



a(9) b(9) b(9) 

a{9) b{9) b{9) 

b{9) a{9) b{9) 

b{9) a{9) b{9) 



(5.15) 
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where a(9) and b(9) are yet to be determined. Requiring that (5.15) is the inverse of 
71/pjT 1 , leads to the following two equations for a{9) and b{9) 



(1 + h(9))a{9) + (n - l)h{9)b{9) = 1 
h{9)a{9) + [{n - l)h{9) + l\b{6) = 



The solution of these simple equations gives 

l + (n- l)h{6) 



a{9) 



and 



KB) 



(5.16) 



(5.17) 



l + n/i(0) ~ K " J l + nh(6) 

Now we can easily obtain the reflection and transmission amplitudes. Its final form is 
easily computed from the formulas above, and we obtain 



R(0) = 



A{9) B{9) A{9) 
B{9) A{9) A{6) 



A{9) 
A{9) 



(5.18) 



A{9) A{9) ... A{9) B{9) 
_A{9) A{6) ... B{9) A{6)_ 

that is, there is a 2 x 2 block-diagonal structure, with diagonal elements A(9) and off- 
diagonal elements B(9), and all remaining elements are equal to A(9). The functions 
A{9) and B(9) are given by 



A{9) 



h{9) 



and 



B{9) = 



l)h{9) 



(5.19) 



l + nh(9) ' v y l + nh(9) 

As one expects, all diagonal elements of R{9) are equal, since there is no reason why 
the reflection of one replica off the defect should be differenct from another. What is 
somewhat surprising is to find out that some of the other elements are equal. This 
means, for example, that a particle of type i coming from the left of the defect has the 
same amplitude of being reflected to the left as a particle of type i, as to being transmitted 
through the defect as any particle different from i. We are interested, of course, in the 
limit n — > 0, where such considerations loose their meaning, though. 



6 Analytical Structure and Correlation Functions 

Since the amplitudes are 27ii-periodic, we can restrict our attention to the physical strip 
—in < Im(#) < in. Let us look at Aiff) initially. As far as < A < 4m/n, there are 
two poles on the physical strip, in the imaginary axis. If A = 4m/n these two poles 
coalesce into one with 9 = —ml, and as A > 4m/ n there are two poles with imaginary 
part — iir/2 and real part different from zero. These correspond to ressonant states, and 
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as A — > oo these poles move to infinity. But since we are interested in the replica limit 
of n — > 0, we see that there are no ressonances for < A < oo. 

We could also analyse the case A < 0, even though A refers to the width of a gaussian 
distribution for the defect disorder. We are looking, therefore, at the model obtained in 
(3.4), as the starting point field theory. In this case we have that for —4m/n < A < 
there are two poles on the imaginary axis. As A decreases and reaches —4m/n, these 
two poles coalesce at in/2. Further decreasing of A gives a real part to these two poles, 
that is, they correspond now to instabilities in the theory. As discussed in [6], one can 
interpret this instability as the emission of pairs of particles from the defect. But in our 
case, since eventually n — > the bound —4m/n — > — oo, that is, the theory becomes well 
defined for all values of the "coupling constant" A. 

We follow the discussion of [6] on how to compute correlation functions. In order 
to do that we rotate the coordinates, t — > ix and x — > —it. Note that since (e,p) = 
(m cosh 9, sinh 9), we have to replcae, accordingly, 9 — > in/2 — 9. We will denote by a hat 
the functions of 9 computed at in/2 — 9. We are interested in computing objects like 

G{x x , h; x 2 , t 2 , . . . ; x n , t n ) = (T[^(x 1 ,t 1 )^(x 2 , t 2 ) . . . $(x n , *„)]) , (6.20) 

where the dynamics of the fields <&(x,t) is governed by the Hamiltonian associated with 
(3.4), which includes the defect interaction term, and T[. . .} denotes the time-ordering 
operator. We can therefore rewrite the correlation function as 

<r[*(*„ (1 )*(* 2 ,t 2 )... *(*„,*„)]> = m^mEi^^^hm . 

(6.21) 

On the left hand side of this equation we compute the expectation value of the string 
of operators using the vacuum of the complete hamiltonian, that is, the bulk part plus 
the defect interaction. On the right hand side the fields are governed by the bulk 
hamiltonian only, and the defect interaction is taken into account by introducing the 
operator V at t — 0. The meaning of this equation using euclidean path integrals (or 
the transfer matrix in the associated statistical mechanics problem) is quite clear, and 
corresponds to the insertion of the perturbing operator V = exp(— / C d ) precisely at 
t = 0. 

Within this framework, we can compute all correlation functions using the form- 
factors of the free bulk theory and matrix elements of the defect operator T>, namely 
(9i,j\V\9 2 ,i), (0\V\9i,i;9 2 ,j), and (9 1 ,i;9 2 ,j\V\0), where \9,i) refers to an asymptotic 
state with rapidity 9 and type i, and so on. Using the path-integral interpretation of 
formula (6.21), we can show that 

(9 1 ,j\V\9 2 ,i) = 2n5(9 1 - 9 2 )f ij {9 1 ) 



7 



(0\V\9 1 ,i;9 2 ,j) = 2it5{9 1 + 9 2 )R ij {9 1 ) 

(9 1 ,i;9 2 ,j\V\0) = 2 7 r5(9 1 + 9 2 )R ij (9 2 ) , (6.22) 



We can compute ($ 2 (x, t)) by using the formulas described in the appendix. In this 
case we have 

(&(x,t)) = ton(& r (x,t) ,) (6.23) 

?1— +0 

where the index r is chosen at will. In the following we denote the coordinates (x, i) by 
z. Using the method described above, we have 

ww; = mm 

= KmY,m\k)(k\V\0) 

k 

= lim — y ^^(0\ ( / ) l\9 1 ,r;9 2 ,r)(9 1 ,r;9 2 ,r\V\0) , (6.24) 

where we have inserted a resolution of the identity, in the asymptotic state basis, be- 
tween (fl and V. The matrix element (O|0^|#i,r; 9 2 ,r) can be computed using the mode 
expansion for the fields 

(O|0 2 |#i v 9 2 r) = 2e~ m< ( coshei+coshe2 ) +imx ( sinhei+sinhe2 ) (6 25) 

and therefore we obtain 

/OO rjfj „ 
— R(9)e- 2mtcoshd 
-oo 2ir 

2mt cosh 9 



A f°° d9 e 

(6.26) 



2m J-oo 2n cosh# 

The fact that this correlation function does not depend on the coordinate x is easily 
understood as a consequence of translation symmetry in a direction paralel to the defect. 
Using that / °° d9 e - 2mtcoshe = K (2mt), and that as u -> oo, K (u) ~ ^J, we find 

-2mt 

(0 2 W) - -fir (6-27) 

as i — > oo. 

As another interesting example, we can compute the correlation function of operators 
located at opposite sides of the defect, that is, (&(zi)Q(z 2 )}, where t± < < t 2 - In this 
case, we need to compute 

/rlf) rift 
^^(0\<f> r ( Zl )\9 1 ,i){9 1 ,i\V\9 2 J)(9 2 J\Mz2M 

h3 

/rift 
^{0\Mzi)\0i,r){e u r\M*)\0)Trr(0i) 
2rr 

2tt coshr ' v ; 



oc 



where t = t 2 — ti, which is always positive. From this expression we can extract the 
asymptotic behaviour of (&(zi)&(z 2 )) as t — > oo. 
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7 Conclusions 



We have studied a simple disordered model that, despite the fact that it is non-integrable 
before disordering, becomes integrable, by the use of the replica trick. The theory thus 
obtained is well defined for all positive values of the width A of the background field 
distribution, and presents no poles as the replica number n — » 0. 

This example has the interesting feature of being a model with a line of defect which 
is amenable by the methods of integrable field theory, with a nontrivial interaction at the 
defect line. 

We used the algebraic method to obtain the reflection and transmission amplitudes 
for finite n, which, together with the matrix elements of the defect operator and the 
form-factors of the free-field theory, were used to compute correlation functions of the 
disordered model. Verifying these results through numerical simulations is an interesting 
problem. 

There are several questions that are worth investigating. We have seen that a given 
non-integrable field theory can become integrable as we disorder it. What about disor- 
dering integrable field theories? It would be interesting to study under what conditions a 
given integrable model will still be integrable after disorder, at least in the relevant limit 
of n -> 0. 

8 Acknowledgments 

We would like to thank L. Moriconi for discussions and encouragement, and A. De Mar- 
tino and G. Mussardo for discussions at an early stage of this work. Financial support 
from CNPq (PROFIX fellowship), CAPES and FAPERJ is gratefully acknowledged. 



9 



Appendix 



In this short appendix we collect some of the formulas used to compute correlation func- 
tions within the framework of the replica trick method. More information and references 
can be found in the review by Bernard [20]. 

We can compute quenched correlation functions from the correlation functions in the 
associated replica model. For example, in order to compute (0(x)) we can show that 



(0{x)) = \ira(O r {x)) . (8.29) 

Similarly we can show that other quenched correlation functions can be computed from 
limits of appropriate correlation functions in the replica model 



(0(x)0(y)) = \im(O r (x)0(y)) , 



(0(x))(0(y)) = \im(O r (x)O s (y)); r^s. (8.30) 



n->0 



These formulas are easily established in the following way. One considers the path-integral 
definition of the quenched correlation functions, for example 

(O(x)O(y)) = J mexp {-^Jh) /P0exp( _g) , (8-31) 

with analogous definitions for other quenched correlation functions (connected too). 
Therefore we can compute the connected two-point correlation function (0(x)0(y)) c 
using the analogous definition of the formula above 

ww) c = / m ex P (—/ *) j^ma '« m 



, 1 5 5 
hm — - 



n->o n 5J(x) 8J(y) 



5J{x) SJ(y) 

,. n / r n \ 



= lim-i^^) + . . . + O n (x)){OM + ... + O n {y))) 

n->0 n 

= \im-(n{O r (x)O r (y)))+n(n-l){O r (x)O s (y)M 

= (O r (x)O r (y)) - (O r (x)O s (y)) r ^ s . (8.32) 

In the above formulas sjfy is the effective action obtained after the integration over the 
random field with n replicas. All the other identities are established in like manner. 
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